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O 

abstract 

> 

' Tensor charge of the nucleon, which will be measured in Drell-Yan processes in polarized 

O ' 

proton-proton collisions at RHIC, are studied in quenched lattice QCD simulation. On the 

VO . 

16 3 x 20 lattice with (3 = 5.7, connected parts of the tensor charge are determined with small 
statistical error, while the disconnected parts are found to be small with relatively large error 



bars. Flavor-singlet tensor charge (5S = 5u + 5d + 5s) is not suppressed as opposed to the 



flavor singlet axial charge (AS = Au + Ad + As) 

43 



X 



1 Introduction 

The parton structure of the nucleon in the twist 2 level is known to be characterized by 
three structure functions gi(x,fi) and hi(x,fi) with x being the Bjorken variable 

and /i being the renormalization scale (see e.g., f\ and gi, which represent the quark- 
momentum distribution and quark-spin distribution respectively, can be measured by the 
deep inelastic lepton-hadron scattering (DIS). On the other hand, hi, which represents the 
quark-transversity distribution, can only be measured in the polarized Drell-Yan processes, 
since it is related to the matrix element of the chiral-odd quark operator. Although such 



1 Invited talks presented at the Australia-Japan workshop on Quarks, Hadrons and Nuclei (Nov. 15-24, 
1995, Adelaide, Australia) and at RIKEN International Symposium on Spin Structure of the Nucleon, (Dec. 
18-19, 1995, RIKEN, Japan) 
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experiments are not yet available, it is planned in RHIC at BNL. Therefore, theoretical 
prediction on hi has some importance. Also, whether there is a "transversity crisis" in 
the first moment of hi as in the case of "spin crisis" in g\ is an interesting question to be 
examined. In this paper, we will concentrate on the first moment of h\{x) and report our 
recent studies using lattice QCD simulations [§]. 

2 Tensor charge versus axial charge 
2.1 definitions 

The first moments of gi(x,fi) and hi(x, /i) are called Ag(/x) and Sq(fi) respectively. They 
are related to the nucleon matrix elements of the axial current and tensor current as follows 

(ps | qi^q \ps) = 2Ms M Aq, (1) 
(ps | qia^jsq \ ps) = 2(s^p v - s v p^) 5q, (2) 

where p^ is the nucleon's four momentum, M is the nucleon's rest mass, and is the 
nucleon's covariant spin-vector. 

In the light cone frame, Aq is interpreted as total quark-helicity in the nucleon, while 5q 
is a total quark-transversity in the nucleon Q: 

A 9(» = / [9i(x,fi) +g 1 (x,fi)}dx (3) 
Jo 

= I [N+(x, n) - N-(x, n) + N+(x, n) - N-(x, y)\dx, 
Jo 

and 

fiqip) — / [hi(x, fi) — h~i(x, fi)]dx (4) 
Jo 

= f [N^(x, n) - N^x, n) - N^(x, fx) + N^x, n)]dx. 

^ 

Here N + (x) (iV_(a;)) denotes the momentum distribution of quarks having the same (oppo- 
site) helicity with the nucleon, while N^(x) (Ni(x)) denotes that having the same (opposite) 
transverse polarization with the nucleon. The quantity with bar is the distribution for anti- 
quark. 
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On the other hand, in the rest frame of the nucleon, Aq (Sq) denotes the quark-spin + 
anti-quark-spin (quark-spin — anti-quark-spin), which can be seen by taking /i = % (// = 
0, v — i) in eq. (0,0). In this frame, estimates by using hadron models are possible. For 
example, relativistic quark models which can reproduce qa = 1-25 correctly give simple 
inequalities; 

| 5u \>\ Au\ , | 5d \>\ Ad | . (5) 

The lower component of the Dirac spinor of the confined quarks plays essential role for the 
above inequalities. 

Drawbacks of such model-calculations are (i) the renormalization scale where the matrix 
elements are evaluated is not clear, and (ii) strange quark contribution, which originates from 
the OZI violating processes, is hard to estimate in a reliable manner. Lattice QCD simula- 
tions can overcome these problems even within the quenched approximation. In particular, 
Am, Ad and As have been studied by two groups 0, || and their results are consistent with 
the recent experimental data on the spin structure of the nucleon |5| . In the quenched level, 
two kinds of diagrams arise: One is the connected amplitude (Fig. la) in which the external 
operator is connected to one of the valence nucleon lines, and another is the disconnected 
amplitude (Fig. lb) where the quark line coming from the external operator is closed by itself. 
The latter gives OZI violating amplitude such as the strangeness content in the nucleon. 
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FIG.l: (a) Connected matrix element of the nucleon. Cross denotes the operator insertion, 
(b) Disconnected matrix element of the nucleon. 



3 Matrix elements on the lattice 

3.1 mass and matrix elements 

We use the standard Wilson's action in our simulation in which two basic parameters read 
(3 = Q/g 2 (g being the bare gauge coupling) and the hopping parameter K. In the following, 
instead of K, we use the "quark-mass" ma = (1/K — l/K c )/2 where a is the lattice spacing 
and K c is the critical hopping parameter at which the pion becomes massless. 

Hadron masses are obtained by the correlation function of composite operators for large 
time t. For example, the nucleon mass is obtained from 

(N(t)N(0)) -> const, x e"" 1 ^, (6) 

with N(t) being the spatially integrated interpolating operator for the nucleon N(t) = 
J cPx{qC^ lr jr :) q)q. On the other hand, matrix element of local operator is obtained as 

R(t) = m) N(o)) ' + { 1 1 } (7) 

Namely the linear slope of R(t — * large) gives matrix element defined on the lattice. 

3.2 some remarks 

In our actual simulation, the following points have been considered [[|. 
(i) To get large overlap of iV(0) with the real nucleon, we use wall source at initial time 
slice t = 0. To do this, we made Coulomb gauge fixing only at t = 0. (ii) To avoid mirror 
source at the final time slice t = tf, we set Dirichlet boundary condition at t/. (iii) To 
compare the matrix element on the lattice with that in the MS"-bar scheme, we multiplied a 
renormalization factor Z(fia) to the obtained data. For the tensor charge, Z calculated for 
fi = 1/a using tadpole improved perturbation reads 

3K 



[1-0.44 a s {l/a)} 



4 



4 Results 



We have done simulations in three cases: 12 3 x 20 {(3 = 5.7), 16 3 x 20 {(3 — 5.7), 16 3 x 20 
(P = 6.0). In this report, we will show the results of the second simulation. We have used 
Fujitsui VPP500 and analyzed 1053 gauge configurations. Each configuration is taken after 
every 1000 sweeps. Three different values of the quark masses K = 0.160,0.164,0.1665 are 
adopted to extract physical quantities in the chiral limit. The statistical errors are estimated 
by the jackknife procedure. 

Hadron masses are extracted by the \ 2 fitting of the data in the interval 5 < t < 10. By 
using the physical hadron masses, m %:P: K = 135,770,498 MeV, one obtains a -1 = 1.42GeV 
(a = 0.14 fm), m = 4.8 MeV and m s = 125 MeV. This also predict the nucleon mass as 
m N = 1.13 GeV. The physical volume of the lattice is V = (2.24) 3 x 2.8 fm 4 . 

X 2 fitting in the interval 6 < t < 11 is also applied for the connected and disconnected 
part of the matrix elements. Shown in Fig. 2 are the data for the correlation function R(t) 
for K = 0.164. For connected u, d contributions (5u con ,, 5d con ) , one can see a clear non- 
vanishing linear slope in 5 < t < 10, while disconnected u — d contribution (5udi S ., 8ddi S .) has 
very small slope evenif it exists. 
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FIG. 2: R(t) as a function of t for medium-heavy quark mass K = 0.164. The black (white) 
circle denotes the connected amplitude 5u con X5d con ), while the black triangle denotes the 
disconnected amplitude 5udi S . = 5d dis . 



6 



Table 1 shows the result of the data extrapolated down to the chiral limit. The result is 
compared with that for Aq in || with the same lattice size and f3. 



Tensor charge Sq (this work) 
16 3 x 20, p = 5.7 
1053 gauge configurations 
(/i 2 = 3 GeV 2 ) 


Axial charge Aq (ref. [3]) 

16 3 x 20, (3 = 5.7 
260 gauge configurations 
(/i 2 = 3 GeV 2 ) 


Su con . = +0.89 (2) 
5u dis . = -0.05 (6) 
5u = +0.84 (6) 


Au con . = +0.76 (4) 
Au d i s . = -0.12 (4) 
Au = +0.64 (5) 


5d con . = -0.18 (1) 
5d dis . = -0.05 (6) 
Sd = -0.23 (6) 


Ad con . = -0.23 (2) 
Ad dis . = -0.12 (4) 
Ad = -0.35 (5) 


5s = -0.05 (10) 


As = -0.11 (3) 


= +0.56 (13) 


AS = +0.18 (10) 



Table 1: Comparison of the tensor and axial charges measured on the lattice. 5u = <5w con . + 
5udis., Sd = 5d con + 5ddi S ., Ss = Ss d i S ., and the same definitions also hold for Ag . 

5 Summary and discussions 

In our simulation, we have found the following. 

1. As for the connected part of the tensor/axial charge, the following inequalities hold: 

| 5u |>| An | , \Sd |<| Ad | . (9) 

This is different from the prediction of relativistic quark models which have universal 
inequality eq.(|). 

2. The disconnected part has still large statistical error and one cannot make definite 
conclusion from our simulation. Nevertheless, there is an indication that (i) the dis- 
connected part is flavor independent i.e., 5udi S . ~ 8d d i S . ~ Ssdis., and (ii) they are small 
but slightly negative. 
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3. Flavor singlet tensor charge = Su + 5d + 5s is not suppressed as opposed to AS: 
5£(3GeV 2 ) = 0.56 (13) <-> A£(3GeV 2 ) = 0.18 (10), (10) 



which implies that there is no "transversity crisis" for the tensor charge. 

Now, what we need to understand is the origin of the smallness of the disconnected part 
of the tensor charge. Since the operator qu^^q (57^75?) is a charge conjugation odd (even) 
operator, 5q(Aq) has a meaning q-spin — g-spin (q-spin + g-spin). This indicates that there is 
a large cancellation of the quark-spin content of the nucleon and the anti-quark-spin content 
of the nucleon. Also, disconnected tensor charge is zero in any order of perturbation theory 
in massless QCD, which might have some relation to its smallness in the non-perturbative 
regime. To clarify the above issue, we are currently collecting more data on 5qdi S . and Aq^js. 
simultaneously with 16 3 x 20 lattice (f3 = 5.7). 
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